Thermal Transport for Many Body Tight-Binding Models 
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We clarify some aspects of the calculation of the thermal transport coefficients. For a tight-binding 
Hamiltonian we discuss the approximate nature of the charge current and the thermal current 
obtained by Peierls substitution which is also identical to the equation of motion technique. We 
address the issue of choosing an appropriate basis for making the Peierls construction for transport 
calculations. We propose a criteria for finding an optimum Wannier basis where the difference 
between the exact current and the approximate one is minimum. Using the equations of motion we 
derive the thermal current for a generalized Hubbard model with density interaction. We identify a 
part which is the contribution from the long range interactions to the heat current. For the Hubbard 
model we derive expressions for the transport coefficients in the limit of infinite dimensions. 



I. INTRODUCTION 

The theoretical description of the thermoelectric re- 
sponse of correlated materials is a fundamental prob- 
lem in condensed matter physics, and a breakthrough ip, 
this area has potential technological useful implications.!!! 
The materials, which have been studied as likely can- 
didates for useful thermoelectric properties, are mostly 
semiconductor alloys and compounds. Materials such as 
Bi2Te3/Sb2Te3 and Si-Ge, which are currently favoured 
for room temperature application, belong to this cate- 
gory. Another class of materials, with potentially useful 
thermoelectric properties, are Ce and La fillecLskutteru- 
dites such as LaFe3CoSbi2 and CeFe3CoSbi2.u Theoret- 
ically these materials have been studied successfully us- 
ing band theoryE Recently Mahan and SofoB have shown 
that the best thermoelectric materials could well be cor- 
related metals and semiconductors (i.e., rare earth inter- 
metallic compounds). The development of the dynami- 
cal mean field theory (DMFT)[for reviews see Refs. 4, 5] 
has allowed new studies of the effects of correlation on 
the thermoelectric response using this method on model 
HamiltoniansBEHa More recent combinations of band the- 
ory and, many- body methods such as the LDA+DMFT 
methodlJfor reviews see Refs. 10, 11] or the LDA++ 
methocO offers the exciting possibility of predicting the 
thermoelectric properties of materials starting from first 
principles JI3 This revival of interest in the thermoelectric 
response motivates us to re-analyze in this paper the fol- 
lowing issues: (1) what is the form of the thermal current 
and the charge current which should be used in realistic 
calculations, and (2) how it should be approximated in a 
DMFT calculation. 

The first question is subtle for two reasons. First, as 
noted early on by Jonson and MahanJlJ the electronic 
part of the thermal current operator contains a quadratic 
and a quartic piece (if the electron-electron interaction is 
non-local) in the electron creation and annihilation oper- 
ators. The contribution of this quartic interaction term 
to the,current has continued to be the subject of discus- 
sion.llj Second, while the form of the thermal current and 
the charge current in the continuum is unambiguous, and 



can be calculated using Noether's theorem,E3c3 DMFT 
calculations require the projection of these currents on 
a restricted lattice model. This involves the computa- 
tion of complicated matrix elements, and in practice an 
approximation which is analogous to the Peierls substi- 
tutionEJ for the electrical current is carried out. It is well 
known that the results of this, construction depend on 
the basis set of orbitals used.EEl This raises the practical 
question of how to optimize the basis of orbitals to be 
used in transport calculations. 

The second question is subtle due to the presence of 
interaction terms in the current. This raises the issue of 
how it should be simplified in the evaluation of the var- 
ious current-current correlation functions and the trans- 
port coefficients. This owestion was first addressed by 
Schweitzer and CzycholEl and by Pruschke and collab- 
orators!^ who stated that within the relaxation time ap- 
proximation, this term can be expressed in terms of a 
time derivative, and the vertex corrections can be ig- 
nored. In the review of Georges et. ala it was stated that 
the results of Pruschke et. at hold beyond the relaxation 
time approximation in the limit of large dimensionality 
when DMFT becomes exact but no detailed proof of this 
statement was presented. 

The following are our main results. (1) In section II 
we address the question of the optimization of the basis 
of localized orbitals for transport calculations, following 
the ideas of Marzari and Vanderbilt.Eil For completeness 
and for pedagogical reasons we discuss in parallel work 
on the charge current, which is simpler and better un- 
derstoodE3 than the thermal current. Our conclusions 
in this context have applications for the computation of 
Born charges in empirical tight-binding models.E3 (2) In 
section III we derive the form of the thermal current to be 
used in tight-binding models, and its dependence on the 
orbitals, using the equation of motion technique intro- 
duced in Ref. 24. Our final expression differs in one term 
from the results of Ref. 15. (3) In section IV we describe 
in detail the diagrammatic analysis of correlation func- 
tions of the current operators. We demonstrate explicitly 
that in the DMFT limit of the transport calculation, the 
vertex corrections (even for those involving the thermal 
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current) can be completely neglected, thereby justifying 
the current practice used in all previous DMFT work. 



A(r), and with coulomb interaction between them. The 
Lagrangian is given by 



II. CHARGE CURRENT 

We consider a system of electrons in a periodic poten- 
tial V(r), in the presence of an external vector potential 



J 



L = 



2m 



d 3 r0 f (V - ieA(r)) 2 i/j - 
1 



-^(r')V(r) 



d 3 rV(r)^ip 



(1) 



r 



Here ^^(r) and tp(r) are the electron field operators with 
usual anticommutation properties. We have ignored the 
spin of the electrons only to simplify the notation. In- 
cluding spin in the following analysis is quite straightfor- 
ward. In field theory, when both high and low energy de- 
grees of freedom are retained, Noether's theorem provides 
a robust procedure to identify the various currentsO The 
theorem associates with every symmetry of the action 
a conserved charge and a corresponding current. The 
charge current is determined by the invariance of the 
action S = J dtL(t), under U(l) gauge transformation 

given by V(r) -► ^(r)e l *( r ) and V f (r) -> V f (r)e~ #(r) - 
The transformation does not produce any variation from 



the interaction term, and the well known expression for 
the charge current is 



j = / d 3 r^(r) (V - ieAM) ip(r). (2) 

m 



The above expression is gauge invariant. The part which 
is proportional to the vector potential gives the diamag- 
netic current. 

In order to facilitate further discussion we will perform 
the standard Noether construction in the Wannier basis. 
In this basis the action (which includes both low and high 
energy degrees of freedom) is 



S 



dt 
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where t^ m 
(n/i\A(r)\mv), 
r'||n 4 ^4,n3/x 3 



(nfj,\H \mv), p',;,„ 
and 

Here Ho 



(nn\p\mv), A££, 



TJH1...IIA 



r - 



p /2m + V(r) is the non- 



interacting part of the Hamiltonian, /i is the band in- 
dex, and R„ defines the lattice positions. W^(r — R„) = 
(r|n/i) form a complete set of orthonormal Wannier func- 
tions. The creation and annihilation operators satisfy the 
anticommutation relation {c^, cj^} = S nm S uu . The gauge 
transformation of the fermionic field operators is equiv- 
alent to the variation 5c% = i J d 3 r4>(r)i(;(r)W*(r - R„) 
and Sc]^ = —i J ci 3 r</)(r)i/' t (r)H / AI (r — R„). Expanding 



4>(r) about the point R„ and keeping only up to V(f> term 
(which is all we need to construct the Noether current) 
we get 

5c£ = z0(R„K + l V<^IO4, 

n t v 

Sclt = -#(R„)c^-iV^^c^LX, (4) 



where = / d 3 rW*(r - R„)(r - R„)W„(r - R m ) are 
the connection coefficients. The matrix L is hermitian, 
i.e., = We note first that the variation from 
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the interaction term is exactly zero. Next, using the op- 
erator identity [rj,Aj(r)] = we find that the variation 
from the term quadratic in A(r) is zero. To get the 



correct diamagnetic part we make use of [ri,pj] — iSij. 
From the invariance of the action we can identify the 
charge current as 



i - »>V*fR — R ^ >r" 4- ie r tM (i^\? v - \? 1 i lv \ r v - — \^ A*r v 

J — lt / J K Ix m "■»)'■„(-„ l-m t IC ^ l„ y-nl^lm Hii L lm> C m / , A nm c n L m 

nm nm i nm 

= ie^c^C<n M |[Wo(A),r]|mi/). (5) 



nm 



Wo (A) = (p - eA) 2 /(2m) + V(r). This is just equation 
(2) expressed in the Wannier basis. The charge-current 
is related to the electronic polarization operatorEj 

Pel = e^4 /i c^(n/x|r|mi/) 

nm 

fXV 

by dP c \/dt = j. The change in polarization AP G i (which 
is a well defined and measurable bulk quantity, rather 
than polarization itself) between an initial and a final 
state of a sample is the integrated current flowing through 
the sample duringj-an adiabatic transformation connect- 
ing the two statesH3 

Theoretical models of the tight-binding type are ef- 
fective low energy models described in tcajas of those 
bands which are close to the Fermi surface .E3 The ques- 
tion, which is non-trivial and which is still debated, is 
what should be the form of the current for such low en- 
ergy models. The low energy Hamiltonian is obtained 
by eliminating or integrating out the degrees of freedom 
corresponding to the high energy bands. This is easily 
formulated in the functional integral language and the 
procedure generates many interaction terms that are not 
present in the original action. In a Hamiltonian formu- 
lation this is equivalent to making a canonical transfor- 
mation-tip decouple the low energy and the high energy 
sectorsE3 That is, given a full many body Hamiltonian Ti, 
we perform unitary transformation U such that UTi.11^ 1 
is diagonal (for a system of interacting particles, in gen- 
eral, this can be done only approximately) , and then con- 
sider only PU7iU~ 1 P ', where P is the operator projecting 
on the low energy bands. To obtain the expression for 
the current in the low energy sector one has to perform 
the same canonical transformation used to transform the 
original Hamiltonian into the effective Hamiltonian on 
the operator representing the current. In other words, 
we first calculate the current (say, J) for the full theory 



(using the symmetry of the full theory), make the same 
unitary transformation and then project the current on 
the low energy sector of interest. The exact low energy 
current is then given by PU3U~ l P. This method of 
calculating the current for the low energy theory is moti- 
vated by renormalization group ideas. But, to implement 
this in practice is usually a formidable task. However, if 
we consider a system of non-interacting electrons (in a 
periodic potential) with a subset M of bands that de- 
fines the low energy subspace, the low energy current is 
obtained by projecting the full current in eqn. (5) on 
the low energy subspace. This is given by PjP, where 
P = |n/x)(n/x| is the projection operator. We note 

that the calculation of the exact current requires knowl- 
edge of the matrix elements of the position operator in. 
addition to that of Ho (the tight-binding parameters). cB 

Sometimes, to avoid calculating the matrix elements 
of the position operator, one makes the approximation 
known as Peierls substitution. There are two types 
of approximations involved with this procedure. First, 
terms involving the connection coefficients are dropped 
out, and one considers an approximate gauge transforma- 
tion given by 8c^ = «^(R„)c^ and fcj/ 1 = — i0(R n )c^. 
Putting the connection coefficients to zero is equivalent 
to the approximation (nii\v\mv) sa R„j„ m i F for the ma- 
trix elements of the position operator, and (n/tt|p|mz/) = 
im{nn\\Ho, r] \mv) w im(H m — R„)i^ for the matrix el- 
ements of the momentum operator. Second, with this ap- 
proximate gauge transformation, the variation from the 
interaction term is non-zero (though, as already noted, 
it is zero for the exact gauge transformation). However, 
contribution to the current from the interaction term is 
neglected. It will be further assumed that the vector po- 
tential is constant, i.e., A^ n = A<5 nm <5 M „. With these 
simplifications the approximate current (jp) is given by 



jp = ie £ ( R ™ - Rt.)C#<, + e2 E ( R ™ - R ») (( R ™ ~ R «) ' A ) C^C ( 6 ) 



nm nm 
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The second term is the approximate diamagnetic contri- 
bution. The usefulness of jp lies in the fact that it can 
be calculated from the tight-binding parameters alone. 

The construction of the Peierls current in terms of 
the atomic orbitals is a priori not obvious for the case 
when there is more than one atom per unit cell. It is 
worthwhile to clarify this issue here. We will denote 
the atomic wavefunctions by |q;tR„), where a is a sym- 
metry index, R n is the lattice position of a unit cell, 
and R T is the position of the atom r within a unit cell. 
It is desirable to define the Bloch basis wavefunctions 
by |ark) = -±= J2 Rn e - 4k '( R "+ R -) |arR„), though the 

phase factor e ~ 4k RT is quite innocuous for the defini- 
tion of the Hamiltonian matrix H(k) aiTi;a2T2 and for the 
subsequent calculation of the energy bands. The ques- 
tion, whether to keep the phase factor or not, is how- 
ever important for the definition of the Peirls current 
jp(k) Ql T i; a 2 T 2 = ^H(k) 0lTi;0jT2 . It is easy to verify 
that, with the above definition of the Bloch basis, one 
gets the same form for the Peierls current if one consid- 
ers a lattice with one atom per unit cell (for which case 
the definition of the Peierls current is unambiguous), and 



compare it with the same lattice with its period doubled 
(and therefore now with two identical atoms per unit 
cell). 

We will examine the behaviour of the exact cur- 
rent and the approximate one under infinitesimal uni- 
tary transformation = 8 nm 5^ u + (where W 
is antihermitian) of the Wannier functions defined by 
\ n ^) ^ J2,niy^mn\ mu ) ■ The variation of a matrix ele- 
ment = ie(n/i|[Wo(A),r]|m^) of the exact current 
is given by 

0)^ - (j)z + E {o&z - w% . (?) 

fc, 7 

This is the usual transformation of matrix elements of 
operators that remain invariant under unitary transfor- 
mation. In fact, the paramagnetic and the diamagnetic 
parts of the operator j are separately invariant. The 
behaviour of jp is however different. The variation of 

(j-P)nm = * e (R"i — Rn)'nm + e _ Rra)((Rm — Rn) ' 

A)t^ m is given by 



Op)z> - (jp)r m + E {(jp)%wz - wjsap)^} + *eE ( R » - R *) - ( R * - r») 

+e 2 E {(Rfc - Rn) ((Rm - Rfc) ■ A) + (R m - R fe ) ((R m - 

fe,7 

-e 2 E {(Rm - Rfe) ((Rfe - Rn) ■ A) + (R fc - R„) ((R m - R„) ■ A)} W^T m - (8) 

fe,7 



The paramagnetic and the diamagnetic parts of jp are 
both basis dependent operators. 

The basis dependence of jp raises the practical ques- 
tion as to what basis one should choose while making the 
Peierls construction. For example, there have been efforts 
to calculate polarization properties, like effective charges 
of semiconductors, using the empirical tight-binding the- 
ory.E3 In this scheme a natural approximation is the "di- 
agonal" ansatz which assumes that the position operator 
is diagonal in the tight-binding basis with expectation 
values equal to the atomic positions. This is equivalent 
to a Peierls substitution, and the polarization calculated 
with this ansatz is related to Peierls current jp. The ef- 
fective charges calculated in this procedure depends on 
the choice of the underlying Wannier basis. In order to 
improve the results one should first make an appropri- 
ate choice of a basis. One possibility is to use the ba- 
sis of the "maximally localized" Wannier functions that 
was introduced by Marzari and VanderbiltEd This is ob- 
tained by minimizing a functional which measures the 
spread of the Wannier functions. Intuitively, it seems 
plausible that the approximation in which the connec- 



tion coefficients are neglected, will work better in a basis 
where the Wannier functions are mace localized. A sec- 
ond possibility, suggested by Millis,E3 is to choose that 
basis in which the charge stiffness calculated using the 
Peierls current will be closest to the one obtained from 
band theory. We note that this criteria is already satis- 
fied by the Bloch basis in which the effective one-electron 
Hamiltonian is diagonal in the band indices. This can be 
seen easily in the following manner. We consider the sce- 
nario of band theory where electrons are in an effective 
periodic potential. Let denote the single particle en- 
ergy levels. It can be shown that the charge stiffness 
is given by D a0 = £) k/i /(e kAl )(d 2 e kjU / (dk a dk ))H Here 
/(e) is the Fermi function and a, (5 denote spatial di- 
rections. The Peierls current constructed in the Bloch 
basis does not have any interband term since the ba- 
sis is already diagonal in the band indices. The para- 
magnetic part of the current is given by (jp) pa ra,a = 

Ek^P ek c/^°) c k c k- Since the paramagnetic part has 
no interband matrix element, it does not contribute to 
the charge stiffness. The diamagnetic part, given by 
(jp)dia,Q = -Y,^p{d 2 e^/{dk a dk l 3))A j3 (^c^, gives a 
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charge stiffness exactly equal to that obtained from band 
theory. It is possible, though, that there are other bases 
which satisfy this criteria. 

In passing we note that if the matrix elements of the 
exact current j are known by some means, say, from first 
principles calculation, then it is possible to define the 
functional 

= ^2 ( n tAi -jp\™>v) ■ (mv\j -jp|n/z) (9) 

nm 



The optimum basis is the one for which the gradient van- 
ishes. The choice of basis will depend on the vector po- 
tential, but the physical quantities calculated in that ba- 
sis will not. In general, this criteria will give a basis which 
is different from that of the "maximally localized" Wan- 
nier functions. The above method of choosing an appro- 
priate basis is not very useful for doing charge transport 
calculations because to define the method one needs to 
know the matrix elements of the exact current, knowing 
which makes the Peierls construction redundant. How- 
ever, one can use this optimization procedure for doing 
thermal transport calculation. As we will see in the next 
section, the matrix elements of the exact thermal cur- 
rent are quite complicated, and a Peierls formulation of 
the thermal current is desirable (in some suitable basis). 
The rationale for our suggestion is that the basis which 



Here p(r) = tp^ (r)ip(r), and U(r) is the two-particle in- 
teraction energy (Coulomb potential, in our case). The 
second term above, which is formally quartic in the field 
operators, is the contribution to energy current from the 
non-local (ia space) interaction. This term was missed 
by Langer,ll3_but noted in a different context by Jonson 
and Mahan.ta More recently, it has been discussed by 
Moreno and Coleman.Ej 

We have discussed in the previous section that for 
an effective low-energy model any current is obtained 
correctly by projecting the current for the full theory 



and choose the basis which minimizes f2, and thereby 
the difference between the exact current and the approx- 
imate one. Using eqns. [7] and [8] we can calculate the 
variation of f2 under infinitesimal unitary transformation. 
The gradient, defined as G^f m — dQ/dW^, is given by 



I 

optimizes the Peierls construction for electric transport 
will be a good basis for doing the Peierls construction for 
thermal transport as well. 



III. THERMAL CURRENT 

In field theory the energy current (which is same as 
the thermal current, except for the latter the single par- 
ticle energies are measured from the chemical potential) 
is determined from the invariance of the action under 
the transformation of time t — ► t — 4>(v,t). This shifts 
the field operators by Sip = ip<p, and Sip^ = ip'(j). From 
the variation of the action defined in equation (1), the 
energy current (js) is given by 



(11) 



(where both high and low energy degrees of freedom are 
present) on the low-energy bands. To implement this for 
the energy current one has to consider variations of the 
Wannier operators (5c- 1 = </>(R;)c^ + V0^Z ■ v L^Cj and 

5cY = 0(ROc] p + V</>£ J> ct I/ L*f under translation of 
time. If we ignore the terms with the connection coeffi- 
cients, we get an approximate current which is equivalent 
to a Peierls substitution. The same approximate current 
can be derived from the low-energy effective Hamiltonian 
using the equations of motionEJ Although we are empha- 
sizing the importance of the exact low-energy current, in 



= (R m - R„) • (nn\[H , (j - $p)]\mv) + ze((R m • A)R m - (R„ • A)R„) • (nn\{H , (j - $p)}\mu) 
- ie^{(R m - R„)(R fe • A) + ((R m - R„) • A)R fe } 

■ - ip\ki){k-f\H Q \mv) + (nfi\H \k'y)(kj\j-jp\rm>)} (10) 



3E 



% I d 3 r 1 j d 3 r 2 (r 2 - r x ) 



U(n - r 2 ) |-0 t (ri) j o(r 2 )^(ri) - V' t (ri)/9(r 2 )^(i - i) + ^ t (ri) / 9(r 2 )7/'(ri)| 
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practice, calculating the exact thermal current is fairly 
complicated. Therefore, we will restrict the derivation to 
that of a Peierls type of energy current for a generalized 
Hubbard model described by the Hamiltonian 

« = E C c fe + E V^n^, (12) 

ij ij 

juv.ffff' 

using the equation of motion technique. Here n^ CT = 
c la c icr- The local energy density (hi) is given by 

hi = ^ E (C c l£4r + ^r c i^c) :t can be sh ° wn that 

i 



2 Z_^i y v ij,crcr' n zcr n j(T' ~ v ji,a' a n jcr' 'Ha J • 
3 



3 

~ 2 / ■ > ji,a' a y jer' if icr jcr' u jct' ia^ia^ja' jff' iff iff Jff jff iff iff J(T y ' V- 1 - ^ 

j 

/i^,ffff' 



where O = i[7Y,0]. The energy current (js) is related 
to the energy density by the continuity equation hi + V ■ 
3E(i) — 0. We define /i(q) = J2 i e _lq ' Ri /ii, and similarly 
jfi(q)- The Fourier transform of the Wannier operators 
are defined by c k(T = e~ Ri c^ a , and similarly for 

c ka-- Here N is the size of the lattice. Comparing with 
the continuity equation we get the energy current 



k 



E v ^Cff' 



c k' ff^k^'Ck'-^o- 



kk' 



where n£ ct = cj, ' 



k^cr c k'+kcr- The nrs ^ two terms 
(the quadratic part) in the above eqn. are contribu- 
tions to the energy current from the electron hopping 
and from the local part of the interactions. The last 
three terms (the quartic part) are additional contribu- 
tions to energy flow iom the long range interactions. 
Moreno and ColemanE_3 have calculated the quartic part 
using Noether's theorem for classical fields, and their re- 



sult is 



k.crcr' 



k,cr' 



n 



jn 



k,ff' 



. We 



2 E V ^ V L^ y n -K 

k, fXV, CT(7 

want to argue that this result is incorrect. We note that 



for classical fields the issue of correct arrangement of op- 
erators is not present. Indeed, if we could commute the 
third operator with the second in each of the last three 
terms of eqn. (14) we would get the result derived in Ref. 
(15). However such commutation will generate an addi- 
tional term ^ ^k^k%r e k"-k4:> c k',ff- Thus, proper 

kk ,^iy,{7 

arrangement of operators is important to get the correct 
form of the energy current, which is naturally captured 
in an equation of motion technique but not while using 
Noether's theorem for classical fields. 

The heat current (jq) is related to the energy current 
by 3q = 3e ML where \x is the chemical potentials The 
chemical potential enters only to shift the single particle 
energies, i.e., right hand side of eqn. (14) gives the heat 

current with the re-definition O = i\H — fjj\f,0], where 
TV is the total particle operator. 



IV. TRANSPORT COEFFICIENTS 

In this section we will examine in detail the derivation 
of the correlation functions of the current operators. We 
will consider only the Peierls type of (charge and ther- 
mal) currents to keep things analytically tractable. In 
Kubo formalism the correlation functions are related to 
the corresponding response functions (the transport co- 
efficients). In the framework of DMFTq it is possible 
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(c) (d) 

FIG. 1: Diagrams in configuration space for thermoelectric 
power. Tii is the interaction term. In (a) and (b) the thermal 
current is a two-point vertex, while in (c) and (d) it is a four- 
point vertex. In the limit of infinite d contribution from (b) 
and (d) can be neglected. 

to derive exact expressions for the transport coefficients. 
The essential simplification in the limit of infinite dimen- 
sions (d) is that the self energy and the vertex terms are 
local. For the single-band Hubbard model, defined by 
the Hamiltonian 

H= (*« C L C .7> + h - c -) + U J2 n i.T n U> 

we will demonstrate that this allows the correlation func- 
tions to be factorized into products of single particle 
Green's functions and their time derivatives. The terms 
that are ignored by such factorization are 0(l/d) smaller 
and can be neglected in the limit of infinite d. Using a 
slightly different approach, the expressions for the trans- 
port coefficients for the Falikov-Kimball model have been 
derived recently.13 

The corporation functions of the current operators are 
defined asE-J 

1 t» 

L ah {iu n ) = —— dre^ T (r r j a (r)j fo (0)), (15) 

piW n V J 

where a, b = (1,2), and ji = j is the charge current and 
J2 = ]q is the heat current. Here V is the volume of the 
system, (3 — 1/ksT is inverse temperature, and iuj n is 
bosonic Matsubara frequency. The transport coefficients 
(that enter the formula for DC conductivity, thermoelec- 
tric power and thermal conductivity) are given by, 

L ab = lim lmL ab (iLu n — > cj + iS). (16) 



For the single band Hubbard model the charge current is 
given by, 

J = e E v k c L<T c k,<r = e 53 i (Rj - R 4 ) Ujc\ a c hai (17) 

k,er (ij) 
a 

and the heat current is given by 

k,<r 

= \ E ( R * - R ;) *« (4, a* - <•*.,'•„) ■ (is) 

(ij) 
a 

Here Vk = Vk£k is the electron velocity. Since the inter- 
action is purely local, there is no contribution from the 
long range interactions. 

The derivation of Lij-pis discussed extensively in 
the literature on DMFT.HE1] In infinite d the-,particle- 
hole vertex becomes momentum independent, Ell and the 
dressed correlation function becomes equal to the bare 
one. This implies the correlation function can be fac- 
torized into a product of single particle Green's func- 
tions, i.e., (T T j(r)j(0)) = -^E M "k^(k,^(k,-T), 

where C/ CT (k, t) = — (7Vck. (T (T)cj i CT (0)) is the fermionic 
Matsubara Green's function. We define the Fourier 
transform Q a (\t,T) = \ J2n e~ 4 "" T <? (T (k, ioj n ), in terms 
of which 

Lu(iuj n ) = 




Q a (k, z) has a possible branch, cut at z — e and SUCk,, z + 
iuj n ) has one at z = e — iuj n !E3 Following MaharO'tj one 
can show 

i^&fMwn + ip n )g a (k,ip n ) = 




—n F (e)A cr (k, e) [<? CT (k, e + iw n ) + £ CT (k, e - w„)] , 



where A a (k, e) = — 2ImG^(k, e) is the spectral function 
and rii?(e) is the Fermi function. After analytic continu- 
ation iuj n — > lo + i5, and after taking the static limit we 
get 

e 2 v-- 2 f°° de ( dn F (e)\ , 2/1 , . n , 
Lu = >vl — ^r L - L )A 2 Jk,e). (19) 

K,CT 

The derivation of L21 is more involved, and is not 
well discussed in the literature. Since the heat cur- 
rent has a part which is a four-point vertex, a priori 
it is not clear whether a factorization of the correlation 
function into products of single particle Green's func- 
tions and their time derivatives is possible. We have 
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Ci,a = -i Y,i Uicim - iU Ci^n^a + i^c lta (and similarly for 
cj a ). We ignore the term with the chemical potential for 
the time being (the result remains unchanged). Due to 
the first term the heat current is a two-point vertex, and 
the corresponding diagrams for L 21 are of the type (a) 
and (b) of Fig. 1. The heat current is a four-point vertex 
due to the second term. The corresponding diagrams are 
of the type (c) and (d) of Fig. 1. In the limit of infinite 
d the scaling of the hopping term is Uj = t*^j\fd (Ref. 

4). This implies that G% - (l/y/d)^ (Ref. 4). One 
can show explicitly that diagrams (a) and (c) are 0{l/d) 
(and higher), and diagrams (b) and (d) are 0(\/d 2 ) (and 



higher). In Fig.l, TLj = U J^- n^n^ is the interaction 
term of the Hubbard Hamiltonian. In the limit of infi- 
nite d the latter drops out, and the factorization of the 
correlation function is possible. In imaginary time 

(TVjQ(r)j(O)) d =°° 

Y d £ V l { < T ^ k , CT (T) C U0)) <rrCk,<r(0)4 i<r (T)) + h.C.} . 
k,<T 

Using £Q{ T ) = (T T £c(T)ct(0)) - S(t) (in imaginary 
time), we get 



L 21 (ius n ) = - 



1 



d) \(3iui n V / 



k,(T I ip n 



iPn + ) ^( k > iPn)Qcr(K ip n + iuj n ) - n Ka 



We drop the second term within braces because it does 
not contribute to ImL 2 i(<^ + iS). The rest is evaluated 
like Ln(iu) n ). It can be shown that 



^-n F (e)A rT (k, e) 



(- + — ) £ CT (k,e + iw n ) 



+ (e y J G„{k,z-iu n ) 



After analytic continuation and taking the static limit we 
get, 



_ e ^ 2 [°° de 



The derivation of L 22 is analogous to that of L 2 i. In 
the limit of infinite rf, (7Vjg(r)jQ) factorizes into prod- 
ucts of (imaginary) time derivatives of single particle 
Green's functions (plus terms which do not contribute 
to Imi 22 (cj)). As in the case of L\\ and L 2 i, the terms 
which are dropped out by such factorization are at least 
0(1/ d) smaller. In other words, 



(TA Q {r)h) ^ 

^ E < { ^k,. (r)cL (0)> (T T c klCT (0)4 iCT (r)) 

k,cr 

-<T T c k , <r (r)4 i<7 (0))<T T c ki(T (0) C j Ci<T (r)) +h.c.} . 



(20) With this simplification it can be shown that 



L 22 (iu n ) = -{\ 



y(3lLO n V , 



k,<j 



The terms in the ellipses do not contribute to ImL 22 (u;). 
Finally we get, 




(21) 

We reiterate the observation made in Ref. (8) that the 
above expressions for the transport coefficients are cor- 
rect for any model with local interaction (for which Eq. 



[18] is correct), in infinite dimensions. 



V. CONCLUSION 

The current (charge or thermal) obtained by Peierls 
substitution or by the equation of motion technique is an 
approximation to the exact low energy current for an ef- 
fective tight-binding Hamiltonian. In particular, the ap- 
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proximate current is not invariant under a unitary trans- 
formation of the Wannier basis. We have suggested a 
simple criteria by which one can choose a set of Wannier 
functions where the difference between the exact and the 
approximate current is minimum. The minimization pro- 
cedure is well defined provided the matrix elements of the 
exact current are known from first principles calculation. 
Using the equations of motion we have derived the ther- 
mal current for a very general tight-binding Hamiltonian, 
correcting the result of a previous work. Finally, using 
the Peierls currents, we have established the correctness 
of known expressions for the transport coefficients for 
the Hubbard model in infinite d. The simplification in 
the limit of large coordination is that the current (charge 



and thermal) correlation functions can be factorized into 
products of single particle Green's functions and their 
time derivatives. These expressions are correct for any 
model with local interaction and in infinite dimensions. 
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